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Introduction 

Over the past several years, transport designers have expressed a renewed interest in 
semi-span model testing, particularly for the design and optimization of high-lift systems. 
Figure 1 shows a schematic of a 2-D high-lift system. This particular geometry consists of 
three elements: a leading edge slat, a main element, and a trailing edge flap. Reference 1 
gives a detailed descreption of the flow around 2-D multi-element airfoils. As shown in F'g. 
1. the wake shed from the leading edge slat travels downstream and subsequently merges 
with the boundary layer developing on the main element to form a confluent boundary layer. 
The confluent boundary layer then travels downstream and merges with the boundary layer 
developing on the flap. The development of confluent boundary layers is highly dependent 
on the Reynolds number. As the Reynolds number increases, the boundary layers developing 
on the slat and main element become thinner. As a result, the location where the slat wake 
and main element boundary layer merge tends to move further aft. Due to such behavior, 
the variation of maximum lift coefficient with Reynolds number has been seen to be highly 
non-linear [2], For this reason, it is necessary to design and test high-lift systems at flight 
Reynolds numbers. 

To meet this challenge, a semi-span testing technique has been proposed for the NASA 
Langley Research Center’s National Transonic Facility (NTF). Semi-span testing has several 
advantages including (?) larger model size, giving increased Reynolds number capability; (??) 
improved model fidelity, allowing ease of flap and slat positioning which ultimately improves 
data quality [3,4] ; and (Hi) reduced construction costs compared with a full-span model. In 
addition, the increased model size inherently allows for increased model strength, reducing 
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aeroelastic effects at the high dynamic pressure levels necessary to simulate flight Reynolds 
numbers. 

The Energy Efficient Transport (EET) full-span model [5] has been modified to become 
the EET semi-span model. The full-span EET model was tested extensively at both NASA 
Langlev Research Center, and NASA Ames Research Center. The available full-span data 
will be useful in validating the semi-span test strategy in the NTF. In spite of the advantages 
discussed above, the use of a semi-span model does introduce additional challenges which 
must be addressed in the testing procedure. To minimize the influence of the sidewall 
boundary layer on the flow over the semi-span model, the model must be off-set from the 
sidewall [3]. Figure 2 shows a schematic diagram of the proposed sidewall mount for the semi- 
span model. The objective is to remove the semi-span model from the sidewall boundary layer 
by use of a stand-off geometry. When this is done however, the symmetry along the centerline 
of the full-span model is lost when the semi-span model is mounted on the wind tunnel 
sidewall. In addition, the large semi-span model will impose a significant pressure loading on 
the sidewall boundary layer, which may cause separation. Even under flow conditions where 
the sidewall boundary layer remains attached, the sidewall boundary layer may adversely 
effect the flow over the semi-span model [6]. Also, the increased model size and sidewall 
mounting requires a modified wall correction strategy [7]. With these issues in mind, the 
semi-span model has been well instrumented with surface pressure taps to obtain data on 
the expected complex flow field in the near wall region. 

This status report summarizes the progress to date on developing the semi-span geom- 
etry definition suitable for generating structured grids for the computational research. In 
addition, the progress on evaluating three state-of-the-art Navier-Stokes codes is presented. 

S emi-span geometry definition 

The definition of the semi-span geometry is needed primarily for developing structured 
grids for the computational research. In addition, the geometry definition can be used for 
developing input to the wall correction codes which will be used in the experimental testing. 
Figure 3 shows the cross-sectional definition of the semi-span model. The wing definition 
was obtained from tabulated theoretical coordinates [5]. The fuselage cross-sections were 
obtained directly from model drawings, which were scanned onto a workstation and dis- 
cretized. These datafiles are stored on the Cray2 (Voyager) supercomputer in the directory: 
milholen/EET. An ascii datafile. README, describes the available stations and their for- 
mat. 


The cross-sections were input into various CAD packages to generate a surface geometry 
definition suitable for use with the GRIDGEN grid generation software. Figure 4 shows an 
oblique view of the resulting surface geometry, including the sidewall stand-off. Several areas 
on the semi-span model had to be approximated, since insufficient information was provided 
from the available sources. These regions included the wing-fuselage fillet, and the wing 
tip. Direct inspection of the semi-span model indicated that the resulting approximations 
provided an accurate representation of the geometry. In addition, the horizontal tail fairing 
present on the semi-span model was added to the database for improved fidelity. Finally, 
several software packages have been utilized to examine the surface slope and curvature 
distributions, to verify the smoothness of the resulting database. 



Computational code evaluation 


The complex geometry of the semi-span model, along with the expected highly three- 
dimensional flow field around the configuration make the flow solver selection even more 
important. With this in mind, three current state-of-the-art Navier-Stokes solvers are being 
evaluated to determine which flow solver would be best suited for the proposed research. 
As discussed above, the semi-span test technique is being implemented in the NTF for the 
design, testing, and optimization of high-lift systems for transport aircraft. The flow fields 
of interest have freestream Mach numbers of approximately .20, which is characteristic of 
take-off or landing conditions. Under such flow conditions, it may seem plausible that the 
flow around the semi-span model is largely incompressible. This however, is generally not 
the case. With high-lift devices deployed, particularly in the leading edge region, localized 
regions of supersonic flow are possible, which may even include shockwaves [1]. 

With these issues in mind, two of the flow solvers chosen for the evaluation solve the 
compressible form of the Navier-Stokes equations. These codes were developed at NASA 
Langley Research Center, and are: TLNS3D [8], and CFL3D [9]. Another logical reason 
for including compressible flow solvers in the evaluation procedure is the advanced state of 
compressible flow solvers. Historically, computational fluid dynamics has been motivated by 
the need to efficiently and accurately obtain solutions for compressible flow problems. As 
a consequence, solutions of the incompressible Navier-Stokes equations have only recently 
become efficient. The third code, INS3D [10], recently developed at NASA Ames Research 
Center solves the incompressible form of the Navier-Stokes equations. All codes have multi- 
block capability, necessary for complex geometries, and several turbulence modeling options 
necessary for the highly three-dimensional flows. 

At low Mach numbers, it may be expected that the performance of a compressible 
flow solver will degrade to the point where either the solver becomes inefficient, or accurate 
solutions may be unobtainable [11]. This inefficiency stems from the fact that the allowable 
local time step, determined from stability analysis, is inversely proportional to the local speed 
of sound [12]. For an incompressible flow, the speed of sound is by definition infinite. As the 
Mach number is lowered for a compressible flow solver, the speed of sound increases, and 
can become quite large. This effectively reduces the allowable time step, which translates 
directly into increased run times. Other problems which can occur at low Mach numbers are 
illustrated with the following example. 

Consider the unsteady compressible continuity equation, written in conservative form: 

| + V.(,v)=0 (1) 

In the absence of high lift devices, the flow around a simple wing geometry is nearly in- 
compressible at low Mach numbers. More precisely, the density is nearly constant, and the 
unsteady term tends to zero throughout the entire computation. This has two direct implica- 
tions. First, the equation set is not guaranteed to be hyperbolic, and thus the time marching 
scheme may become inappropriate. Secondly, many codes use this unsteady term to monitor 
convergence to steady state. If this term is essentially zero, it may be necessary to monitor 
another term, or a global parameter such as lift coefficient. In addition, if the flow field is 
essentially incompressible, it is unnecessary to solve the energy equation simultaneoulsy with 
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the continuity and momenta equations. For this reason, the energy equation can decouple 
from the solution procedure, introducing oscillations, and degrading the accuracy of the code 
[13]. This underscores the need to evaluate the performance of the compressible flow codes 
at low freestream Mach numbers. 

For the evaluation procedure, an experimental data set was chosen which had both a 
simple geometry, and a freestream Mach number similar to that proposed for the EET semi- 
span test technique. Figure 5 shows a partial view of a C-0 grid employed for the first wing 
geometry. The unswept rectangular wing has an aspect ratio of six, and is constructed of 
the NACA 0012 airfoil section. In comparison, the full-span EET model has an aspect ratio 
of 10. The wind tunnel model had a rounded wing tip, which has been modeled in the grid 
shown in Fig. 5. Reference 14 gives details of the available experimental data. The second 
wing geometry is the same wing swept aft 30°. This leading edge sweep angle was chosen 
because it is nearly identical to that of the EET semi-span model. It should be noted that 
a wing tip extension was added to the swept wing so that the wing tip would be parallel to 
the approaching freestream. 

Figure 6 shows a comparison of the computational results obtained from both compress- 
ible flow solvers to experimental data for the following conditions: = .14, a = 6.75°, 

Re = 3.3 x 1 0 6 , and A = 0°. Both computations agree quite well with the experimental data. 
At the inboard station, both computational results are nearly identical. At the outboard 
station the computations show some differences, however, the predicted pressure distribu- 
tions are still in good agreement with the data. For brevity, other stations across the wing 
will not be examined, but similar agreement with data was obtained. This overall agreement 
over the entire wing surface is illustrated by comparing the computed spanwise load distri- 
butions to the experimental data, as shown in Fig. 7. Here tj represents the non-dimensional 
semi-span fraction, c n the sectional normal force coefficient, and Cat the total wing normal 
force coefficient. Both computations agree quite well with the experimental data. 

Figure 8 compares the computed pressure distributions to experimental data at two 
spanwise locations for the swept wing geometry. The conditions for this case are: 

Mqo = -14, a = 6.75°. Re = 3.3 x 10 6 , and A = 30°. At both stations, the computations 
show good agreement with the data. The leading edge suction peaks have been consistently 
underpredicted across the wing surface. This may raise concern about the capability of the 
compressible solvers to predict the flow around swept wings at low Mach numbers. It should 
be noted however, that the reflection plane used during the wind tunnel experiment was 
small in comparison to the wing [14]. In fact, for this particular sweep angle, the tip of the 
wing was well behind the reflection plane. Given this situation, symmetry at the root plane 
of the wing would be unlikely [15]. 

These two test cases have shown that both compressible flow solvers are capable of 
accurately predicting the low speed flow around the 3-D wing geometries. The obtained 
solutions showed no signs of oscillations, indicating that no obvious numerical difficulties 
were introduced by t he low freestream Mach number. 

Due to difficulties with boundary condition implementation. INS3D could not be run for 
the C-0 grid topologies used above. To demonstrate the codes capability, two-dimensional 
computations were performed for a NACA 4412 airfoil. For comparison, CFL3D was run 
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in a 2-D mode, with a freestream Mach number of .20. Figure 9 compares both computed 
solutions to experimental data [16] . for the following conditions: a = 6.10°, Re = 1.8 x 10 b . 
Both computations qualitatively agree well with the data. The INS3D solution shows better 
agreement than the CFL3D solution, which may be due to the artificially high freestream 
Mach number for the CFL3D solution. However, the INS3D solution time was approximately 
four times longer than the CFL3D solution. This is attributed to the fact that INS3D does 
not presently have a multigrid capability. 


Future work 


INS3D computations 

At present, a new version of INS3D has been obtained. The capability of the code 
to compute 3-D solutions using C-0 grid topologies will be evaluated. Once this has been 
completed, a computational code will be chosen which will be used for the computational 
research. 

Flat plate boundary layer calculations 

A flat plate boundary layer code [17] will be used in conjunction with experimentally 
measured sidewall boundary layer profiles to determine the virtual origin of the sidewall 
boundary layer. This will provide velocity profiles which will be imposed at the far-field 
upstream boundary of the computational domain to simulate the sidewall boundary layer. 
The flate plate computations will also be useful in examining what influence the wind tunnel 
sidewall temperature has on the development of the sidewall boundary layer. In addition, real 
gas effects will be examined, to determine how the developement of the sidewall boundary 
layer is effected. 

Semi-span grid generation 

The geometry database discussed above is currently being used with GRIDGEN to 
generate multiblock structured grids for the computational research. 

Free-air computations 

The first phase of computations for the semi-span model will focus on the prediction of 
the “free-air” [6] flow field around the model. This will serve two distinct purposes. First, 
the computations will provide a baseline which will be used to determine what influence the 
sidewall boundary layer has on the flow field. The free-air computations will be validated by 
making direct comparison with the available experimental data. Comparisons will be made 
to the wing pressure distributions and the available force and moment coefficient data. Once 
validated, the code would be used as a diagnostic tool to examine the details of the solution. 

The second purpose of the free-air calculations is to provide detailed information which 
will not be easily obtainable in the wind tunnel testing. Such detailed information is not 
presently available for this geometry, and would be quite useful in the upcoming experi- 
mental investigation. This data will also provide the researchers an independent method 
for evaluating the quality of their data, particularly relating to the wall interference and 
correction strategy. 
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Viscous sidewall modeling 

The viscous sidewall computations are the primary focus of the research project. The 
goal is to examine in detail, what influence the sidewall boundary layer has on the flow over 
the semi-span model. This determination will be made by making direct comparison with the 
free-air computations, and available experimental data. The comparisons will be made on 
several levels. The most general comparisons will be between force and moment coefficients, 
and will give an overall assessment of the influence of the sidewall boundary layer. Detailed 
comparisons between surface pressure distributions, over the entire model, will give insight 
as to how the flow field has been modified by the presence of the sidewall boundary layer. 
The finest level of comparison would examine details of the wing and fuselage boundary 
layers, and how they have been influenced by the sidewall boundary layer. 

The predicted behavior of the sidewall boundary layer is of particular interest. Due 
to the difficulty in obtaining off-surface measurements in the cryogenic environment, little 
experimental data can be obtained to characterize the behavior of the sidewall boundary 
layer. In this particular aspect, the computations will give insight into the inherent juncture 
flow present in semi-span testing. If the predicted characteristics of the sidewall boundary 
layer are not acceptable, such as the formation of horseshoe vortices, techniques to improve 
the flow quality can be computationally examined (such as suction, tangential blowing, or 
vortex generators). 
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Fig. 1: Schematic of ‘2-D high lift system. 
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Fig. 3: Cross-sectional geometry definition of EET semi-span model. 
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Fig. 4: Surface geometry definition of EET semi-span model. 
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Fig. 5: Partial view of C-0 grid topology. 
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Fig. 7: Comparison of computed spanwise load distributions with experimental data 
t M„ = .14. q = 6.75°. Re = 3.3 x 10 6 , A = 0°) . 




Fig. 8: Comparisoi 
(Moo = -14, 
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Fig. 9: Comparison of computed pressure distributions with experimental data 
(NACA 4412, a = 6.40°, Re = 1.8 x 10 6 ) . 


